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ABSTRACT 

We predict the metallicity probability distribution function (PDF) of long gamma-ray 
burst (GRB) host galaxies at low-redshifts (z < 0.3) when GRBs occur only in low- 
metallicity environment, assuming empirical formulations of galaxy properties. We 
discuss contribution of high-metallicity galaxies to the cosmic rate of low-metallicity 
GRBs, taking internal dispersion of metallicity within each galaxy into account. As- 
suming GRBs trace low-metallicity star formation < Z CT n\ 12+logio(0/H) = 8.2, we 
find that GRB host galaxies may have systematically higher-metallicity than that 
of GRB progenitors. Furthermore, we expect > 10% of the host galaxies to have 
12+logio(0/H) > 8.8, if galaxies have internal dispersion of metallicity comparable to 
that observed in the Milky Way. Our results show that the low-metallicity scenario of 
GRB progenitors can be reconciled with the recent discoveries of the high-metallicity 
host galaxies of GRBs. We also show possible bimodality in the host metallicity PDF 
that results from the single progenitor model of GRBs. If found in future observation, 
the bimodality can be a clue to constrain the nature of GRB progenitors. 

Key words: gamma-ray burst: general - galaxies: abundances. 



1 INTRODUCTION 

It is now broadly agreed that some of long gamma-ray bursts 
(GRBs) originate in core-collapse supernovae (CC SNe). 
However, not all CC SNe produce GRBs, and the criteria 
that discriminates GRBs from general CC SN is one of the 
most outstanding questions about GRBs. 

Some theoretical studies on the origin of GRBs using 
stellar evolution models suggest that a low-metallicity may 
be a necessar y condition for a GRB to occur (Z < a few 
x O.lZn, e.g-lMacFadven fc Wooslevll 19991: lYoon fc Langeil 
120051 : lYoon et al.ll200tj : IWooslev fc Hegerll200d ). It has also 
been suggested from observations that metallicity distribu- 
tion of GRB host galaxies at redshift z < 0.25 is signifi- 
cantly biased towards low metallicities compared to the ex- 
ectation when GRBs are unbiased tracers of star formation 



Stanek et al.ll2006l : lModiaz et al.ll200Sl ). 

Furthermore, some observations suggest that the GRB 
host galaxies are systematically fai nter and smalle r than 
those of the core- collapse SNe ihe Floc'h et all 120031 : 
Fruc hter. et al ] l200d : ISvensson et alJl201Cl ). indicating that 
the GRBs may preferentially occur in low metallicity en- 
vironment, because fainter and lower mass galaxies gener- 
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ally have lower metallicities. These interpretations have also 
been supported by other theoretical studies using the mod- 
els of galaxies (e.g. IWolf fc Podsiadlowskil [20071 : iNiino et all 

Eon]). 

However, recent discoveries of high-metallicity host 
galaxies of some GRBs cast doub t on the requirement of 
low-metallicity i n GRB occurrence jLevesaue et al.ll2010al lbl: 
lHan et al.ll2010l : lHashimoto. et ajj|2010| ). On the other hand, 
it should be kept in mind that GRB host galaxies may 
have different metallicity from that of GRB progenitors, be- 
cause a galaxy is not a chemically homogeneous object. To 
decide whether the discoveries of the high-metallicity host 
galaxies are consistent to the low-metallicity requirement or 
not, we need to quantitatively consider contribution of high- 
metallicity galaxies to the cosmic rate of low-metallicity star 
formation. 

In this study, we predict the metallicity probability dis- 
tribution function (PDF) of GRB host galaxies at z < 0.3, 
assuming the low-metallicity requirement and taking the 
metallicity dispersion within each galaxy into account. We 
discuss whether the expected rate of low-metallicity GRBs 
in high-metallicity galaxies is significant to explain the ob- 
servations or not. We use empirical formulations of galaxy 
properties, and assume GRBs trace low-metallicity star for- 
mation < Z crit : 12+logio(0/H)= 8.2. In § [2] we describe 
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empirical formulations of galaxy properties and GRB rate 
models which we use in this study. In § [3] we show the 
expected metallicity and mass distributions of GRB host 
galaxies and discuss their implications. In §[4] we also discuss 
how the expected metallicity distribution changes if there is 
a correlation between metallic ity and star formation rate 
of a galaxy, as it is claimed in iMannucci et~ai1 j20ld ). We 
summarize our conclusion in § [5] 



2 MODELS 

2.1 Galaxy Properties 

To compute the expected metallicity PDF of GRB host 
galaxies at z < 0.3, we assume empirical formulation of stel- 
lar mass function, mass-star formation rate (Af*-SFR) re- 
lation of galaxies, and mass-metallicity (Af+-Z ga i) relation 
of low-redshi f t gala xies. Our approach is similar to that of 
IStanek et al.l |2006), who have calculated expected metallic- 
ity PDF of GRB host galaxies when GRBs trace star forma- 
tion without any metallicity dependence. However, we step 
further to include the low- metallicity preference of GRBs 
considering the chemical inhomogeneity within galaxies. 

It should be noted that A/* is calibrated using dif- 
ferent initial mass functions (IMFs) in different studies. 
In this study, we assume conversion among stellar mass 
scales for different IMFs as: Ms a ip Ctcr = 1.43 x Afdictsaipctcr 

= 1-5 X MKro u pa = 1-8 X Mchabrie r = 1-8 X A/ B G03 

JSalpeterl [l955|; iBell fc de Jond l200ll; iKroupa fe Burkertl 
l200ll ; IChabrieri(2003l ; lBaldrv fc Glazebrookll2003l 'l. Hereafter 
stellar masses are scaled for Salpeter IMF, unless otherwise 
stated. 

We u se empiric al form ulations of the stellar mass 
function (IBell et all |2003|) . and the A/*-SFR rela- 
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a fit to the observation by 
of low-redshift (z < 0.3) late- 
type galaxies. We assume the di spersion of the A/^ -SFR 
relation to be ~ 0.3 dex following IStanek et al.l (|2006l ). Us- 
ing the mass function and the A/^-SFR relation, we com- 
pute cosmic SFR density as a function of M*: psfr(A4*) 
[MQyr -1 Mpc -3 dex -1 ]. These models are shown in Fig. [1] 
We only consider galaxies with logio Af* > 8.0, which corre- 
sponds to the lowest-mass of GRB host galaxies ever known. 
Both of the stellar mass function and the A/*-SFR relation 
may suffer from selection effects of the galaxy sample. To 
demonstrate the uncertainty, we also use psFiifAf*) derived 
from observation of galaxies of all-type l|Drorv fc Alvare3 
120081 . DA08). 

Various methods have been proposed to measure gas 
metallicity of gal axies, but they do not always agree with 
each other (e.g., iKennicutt et alj 120031 ; iKewlev fc Ellison] 
2008) . When we discuss metallicity of galaxies, we use metal- 
licities calibrated with iKobulnickv fc Kewlevl ([200J, here- 
after KK04) method in this study. We use the A/*— Zg a i 
relation with t he K K04 calibration presented in Eq. 8 of 
ISavaglio et al.l (|2005l . hereafter S05 relatio n). We assume 
the d ispersion of the relation to be 0.1 dex dTremonti et al.l 
l2004h . 

When we consider the dispersions of A/^-SFR relation 
and Af*-Z ga i relation, we assume that the offset from A/*- 
Zgai relation and that from Af*-SFR relation are indepen- 
dent of each other. This assumption is supported by the 
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Figure 1. Empirical formulations of stellar mass fun ction, M*— 
SFR relation, and SFR density for late-type galaxies (Bell et all 
1 20031 ; [Brinchmann ct al. 2 0041) , Those for all types of galaxies 
(DA08) arc plotted together. In panel (b), the typical dispersion 
of galaxy SFR around the relation is shown. 



obser ved no correlation be tween Ha equivalent width and 
SFR l|Tremonti et al.ll2004T l. However, it is also clai med that 
those offsets are correlated (|Mannucci et al.ll2010l ). and we 
discuss the case where of the offsets are correlated in § [4] 



2.2 GRB Rate and Internal Dispersion of 
Metallicity within Each Galaxy 

Observations of nearby galaxies, including the Milky Way 
(MW) and the Magellanic clouds, show that the galaxies 
have internal dispersion of metallicity within them (~ 1 dex 
in M W and ~ 0.3 dex in the Larg e Magellanic Cloud (LMC), 
e.g- lRolleston et alJlibocl . I2002T ). Furthermore, there is a ~ 
0.4 dex variation of 12+logio (O /H) among HII regions in the 
host galaxy of GRB 980425/SN 1998bw w hich is comparable 
to 3a error of the metallicity calibration (|Christensen et all 
2008). To demonstrate effects of the chemical inhomogene- 
ity, we assume metallicity of SFR in a galaxy has a log- 
normal distribution with dispersion az,int around Z ga i, al- 
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Figure 2. The log-normal models of the internal metallicity 
dispersion within a galaxy is shown (oz,\nt = 0.1,0.3 and 0.5, 
solid curves), together with the metallici ty distributions of the 
HII regions and the young B-type stars l|Afflerbach et al.lll997l : 
iRoIIeston et al.| [2000, red dashed and bl ue dotted, respectively) . 
We note that the metallicit ies measured in lAfflerbach et al.l l|l997f l 
and lRolleston et al.l (2000) may be inconsistent with KK04 metal- 
licities. 



though metallicity distribution of star forming gas within a 
galaxy is hardly understood. 

In Fig. [21 we plot the log-normal models with Z ga i = 8.9 
and az.int = 0.1,0.3 and 0.5. For comparison, we plot 
the observed metallicity distributi on of HII regions and 
young B-type stars in the MW (|Afflerbach et all 1 19971 ; 
IRoIIeston etai]|2000h together, although they do not nec- 
essarily represent overall metallicity PDF of star forming 
gas in the MW. Both of the HII regions and the B-type 
stars typically have 12+logio (O/H) ~ 8.9, while ~ 10% have 
12+logio(0/H) < 8.2, which is roughly comparable to the 
case of az,int = 0.5. Although the models and the obser- 
vation do not agree with each other in the high-metallicity 
end, we note that the high-metallicity end of the internal 
metallicity distribution has small effect on our results. 

In this study, we assume GRB rate is proportional to 
low- metallicity SFR with Z < Z C nt, and express /?grb oc 
£GRB(Z g ai)xSFR, where e G RB(Z g ai) = SFR z <z crlt /SFR is 
a GRB efficiency in a galaxy. We do not need to determine 
the proportionality constant, because our interest is in the 
metallicity PDF of GRB host galaxies. Under the assump- 
tion of the log-normal metallicity distribution, the GRB ef- 
ficiency can be written as: 



£GRB(Z g ai) = 0.5 x erfc[ 



-log 10 (Z cri t/Zgal) j 
y2<7z,int 



(1) 



Here we assume Z cr it to be 12+logi (O/H) = 8.2, which 
roughly corresponds to 0.2-0.4Zq . The expected metal- 
licity distribution of GRB host galaxies is />grb (Z ga i) oc 
eGRB(Zgai)/3sFR(Z g ai), where /3sFR(Z ga i) is the psfr(M*) 

projected to Z ga i axis using th e M*-Z g ai relati on. 

It should be noted that iKocevski et al. (2009) per- 
formed similar investigation to that in this study to calcu- 
late the mass PDF of GRB host galaxies when GRBs trace 
low- metallicity SFR, and claimed Z cr it > O.5Z0 which is 
contrary to our assumption. However, they assumed that 
Rgrb = when Z ga i > Zcrit without taking the internal dis- 



persion of galaxies into account, and their conclusion might 
be affected by this assumption. 



3 THE METALLICITY PDF OF GRB HOST 
GALAXIES 

The predicted metallicity PDF of GRB host galaxies is 
shown in the left panels of Fig. [3] The model without metal 
cutoff (i.e. £grb (Z g a i) = 1.0) is consistent to the results of 
IStanek et al.l f200fj), and it shows that more than 50% of 
low-redshift star formation takes place in high-metallicity 
galaxies with 12+login (O/H) > 8.8, which is much higher 
fraction than high-metallicity galaxies in observed GRB host 
galaxies. 

Now we consider the effect of the metal cutoff on the 
metallicity distribution of the host galaxies. The contribu- 
tion of Z ga i > Zcrit galaxies is not zero due to the effect 
of the internal dispersion. The results with az,int = 0.1,0.3 
and 0.5 are shown in Fig. [3] In the cases of az,mt J? 0.3 dex, 
more than 50% of GRB host galaxies have Z ga i > Zcrit j sug- 
gesting that the progenitor metallicity can be systematically 
different from the host metallicity. 

The contribution of the high-metallicity galaxies 
(12+logio(0/H) > 8.8) to the cosmic GRB rate is equivalent 
to that of Zgai < Zcrit galaxies when they have <Jz;mt = 0.5 
dex. In the case of oz,iv& = 0.5 (0.3) dex, roughly 25% (5%) 
of the host galaxies have 12+logio (O/H) > 8.8, suggesting 
the hypothesis that GRBs occur only in low-metallicity en- 
vironment does not contradict to the recent observations of 
high-metallicity host galaxies of GRBs. We note that the 
prediction of the DA08 model is not largely different from 
that of the late- type galaxy model. The expected mass PDF 
of GRB host galaxies is shown in the right panels of Fig. [3] 
in the similar manner to the metallicity PDFs. 

We have assumed that az,int is same among all galaxies 
in discussions above. However, the internal metallicity dis- 
persion within galaxies is not well known, and it is possible 
that galaxies with different M* (or Z ga i) typically have dif- 
ferent <7z,mt- In fact, the LMC has smaller internal disper- 
sion o f metallicity than that in MW (e.g., IRoIIeston et al.l 
|2002| ; ICionil 12009 1 In that case, the expected metallicity 
distribution would be different from those discussed above. 
However, the relative contribution of 12+login (O/H) > 8.8 
galaxies compared to the Z < Z cr it galaxies would be similar 
to the case of o"z,int 0.3 dex when high-metallicity galaxies 
have az,int compareble to that observed in the MW, regard- 
less of az,int in 12+login(0/H) ^ 8.8 galaxies. 

To demonstrate a case in which crz.int correlates with 
M^, we employ a toy model of dependent az.int, 



o-z,i„t(M*) = 0.2 x (logM* - 8.0). 



(2) 



The results are shown in Fig. [3] It is interesting that the Z ga i 
PDF has multi-peak distribution in this model, although we 
consider only one population of GRB progenitors. 

This bimodality can be explained as follows. In the case 
of crz,int = 0.5 dex, pgrb (Z g ai) is approximately constant be- 
tween 12+logio(0/H) = 8.2 and 9.0 (see left bottom panel 
of Fig. O. If (1) high-metallicity galaxies (12+logi (O/H) 
> 9.0) typically have az,int ~ 0.5 dex and (2) az,int is pos- 
itively correlated with Z ga i, galaxies with 12+login(0/H) 
< 8.2 or > 9.0 have similar pgrb (Z ga i) to that in the 
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Figure 3. The Z ga i PDF {left panels) and M+ PDF {right panels) of GRB host galaxies predicted in our model. The top panels show 
the normali zed PDF, while the bo ttom panels show the -Rgrb weighted metallicity function of galaxies. Observed GRB host galaxies 
at z < 0.3 (Lcvcsquc ct al. 201Q3, L10) are shown together. The PDFs predicted without the effect of Z CI [ t are plotted with red solid 
lines, while results with crz.int = 0.5, 0.3 and 0.1 dex are shown with dashed (magenta), dotted (blue), and dot-dashed (dark blue) 
lines, respectively. The green double-dotted dashed line represents the results with DA08 models. The results in the case where az,\nt 
correlates with M* are shown with double-dotted line (brown). 



case of CTz.int = 0.5 dex, while smaller cgrb (-^gai ) in 8.2 < 
12+logi (O/H) < 9.0 produce a dip in the metallicity PDF 
between 12+logio(0/H) = 8.2 and 9.0. Note that e G RB and 
o~z. int are positively correlated when Z ga i > Z CI i t . Thus it 
is possible that bimodal distribution of Z ga i appears from 
a single population of GRB progenitors, when the condi- 
tions (1) and (2) are fulfilled. We should keep in mind that 
bimodality of GRB host galaxy population does not neces- 
sarily mean bimodal nature of GRB progenitors. 



4 CORRELATION BEWEEN Z gai AND SFR 

In the previous sections, the dispersions of the M*-Z ga i re- 
lation and the M*-SFR relation are treated independently. 
However, it is recently claimed that galaxies with higher- 
SFR tend to have lo wer- Z ga i compared to l ower-SFR galax- 
ies wi th si milar M+ jMannucci et alj|20ich . iMannucci et al.l 
|201ll ) and iKocevski fc WestJ (|2oTof ) have investigated the 
effect of the SFR-Z ga i correlation on the M*-Z ga i relation 
of GRB host galaxies. 

In this section, we use M*-SFR-Z ga i relation 
l|Mannucci et al.l l201ll . hereafter Mil relation) instead of 
S05 relation, to investigate how the metallicity PDF of 
the host galaxies is changed if there is a SFR-Z ga i cor- 
relation. It should be noted that Mil relatio n studied in 
IMannucci etafl (|20ld . l201ll ) is calibrated with Naga o et al.l 
( 20061 . hereafter N06) method. Hence we can not directly 
compare predictions of Mil relation with predictions and/or 
observations with KK04 calibration. We project the mass- 
SFR relation described in ^ 12. 1 I to M*— Z &a i plane using Mil 



relation, and compare it with S05 relation. In Fig. [4] one 
sees discrepancy between the two relation. 

To make consistent comparison between the two mod- 
els, we assume ad hoc conversion between the two metallicity 
calibration: 



QKK04 = 



ON06 + 0.25(O!N06 
CVN06 — 0.33(aN06 



8.0) +0.1 (q N 06 < 8.4) 
8.4) + 0.2 (a N06 ^ 8.4) 



(3) 



where a represents 12+logio(0/H). Converted with Eq. O 
the projected M*-SFR relation agrees with S05 relation in 
0.04 dex (Fig.g). 

The metallicity and mass PDFs predicted using Mil re- 
lation are shown in Fig. [5] Although the host galaxies have 
lower- metallicity by < 0.1 dex compared to the case of S05 
relation depending on oz,mt, Mil relation alone does not 
make the matellicity PDF consistent to the current sample 
of GRB hosts without further metallicity effect. The metal- 
licity PDFs with the metallicity cut-off may agree with the 
observations, as well as in the case of S05 relation. 

The predicted host galaxies in the case of az,mt = 0.5 
have smaller M* compared to that for S05 relation, while 
the host galaxies in the case of <7z,int =0.1 have larger M*. 
As a result, the mass PDF is less sensitive to the change 
of crz.int. This is because SFR and Z ga i correlate stronger 
when Mi, is smaller in Mil relation. Once a galaxy sample 
is weighted with SFR, Mil relation makes M*-Z ga i relation 
steeper, and hence a difference in Z ga i corresponds to smaller 
difference in M* with Mil relation compared to the case 
with S05 relation. However, we note that the predicted mass 
PDF is strongly dependent on the low-metallicity tail of Mil 
relation, which is still highly uncertain. 
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Figure 5. Same as the top panels of Fig. [3] but Mil relation is assumed instead of S05 relation. The results with S05 relation, which 
have been shown in Fig. [3] are plotted together (grey double-dotted dashed lines). 




Figure 4. The M-^—Z^i relat i on of low-redshift galaxies [z < 
0.3) presented in lSavaelio et all ||2005| . S05 rel ation, the red solid 
line) is compared with the mass-SF R relation jBrin chman n et al.l 
2004) projected by Mil relation jMannucci et al.l l201ll . blue 
dashed line). The conversion of the projected M+— SFR relation 
with Eq. [3] is shown with a cyan dotted line. 



5 DISCUSSIONS 

We have predicted the metallicity and mass PDFs of GRB 
host galaxies, assuming empirical formulations of galaxy 
properties and the model of GRB rate in which GRBs oc- 
cur only from low-metallicity stars (< Z C iit- 12+logio (O/H) 
= 8.2). Our results show that > 50% of GRB host galax- 
ies can have z? ga i > Z cr it, and high- metallicity galaxies 
(12+logio(0/H) > 8.8) may have significant contribution to 
cosmic GRB rate. This means that metallicities of GRB host 
galaxies may be systematically different from those of GRB 
progenitors, and the low-metallicity scenario can be recon- 
ciled with the observations of high-metallicity host galaxies 
of GRBs. 

For some GRBs, met allicities are measu red at the po- 
sitions of the bursts (e.g.. iModiaz et al1l200ct ). and the host 
galaxy of G RB 020819 has high-me tallicity at the position 
of the burst l|Levesaue et al ] |2010al ). However, it should be 
noted that the positioning error of GRB 020819 is roughly 



5 kpc (jjakobsson et al.l [2005 s ) . and there might be chemi- 
cal inhomogeneity in the error circle. We need more precise 
localization of GRBs to draw robust conclusions, although 
it is also difficult to specify what precision is required. The 
required precision is dependent on the mixing process of 
inter-stellar medium which is not well understood. 



Although we have formulated £GHB(^ g ai) motivated by 
the probable chemical inhomogeneity within each GRB host 
galaxy, similar formulation of eGRB(Zgai) may be obtained 
considering other effects (e.g. moderate low-metallicity pref- 
erence of GRB occurrence without sharp metallicity cutoff). 
It is currently difficult to distinguish what effect constructs 

eGRB(Zgal)- 



We have shown that multi-peak distribution of the 
metallicity of GRB host galaxies can be produced by a single 
population of GRB progenitors, when o"z,mt positively cor- 
relates with M*. If observed, the bimodality can be a clue 
to investigate the nature of GRB progenitors. If egrb > 
in high-metallicity galaxies is caused by the nature of GRB 
progenitors rather than properties of galaxies, there would 
be no effect of the o"z,int-Af* correlation. 



Although some results shown in this paper suffer from 
uncertainties about the properties of low-redshift galaxies, 
some important features of the results which we have dis- 
cussed are not dependent on the detail of the modelings (see 
§ |3j . However, we need to understand the actual metallicity 
distribution within young star forming galaxies to make reli- 
able prediction of the exact metallicity PDF. More detailed 
study of the internal structure of galaxies requires different 
approach from that in this study, such as high-resolution 
hydrodynamic simulation and/or spatially resolved spectro- 
scopic observation of large sample of galaxies, and we ad- 
dress this issue to future studies. Future development of our 
knowledge about galaxy properties would provide us with 
more robust predictions about GRB progenitors and their 
host galaxies. 
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